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Summary 

This report describes the ceramics analysis and reliability 
evaluation of structures (CARES) computer program. The 
primary function of the code is to calculate the fast-fracture 
reliability or failure probability of macroscopically isotropic 
ceramic components. These components may be subjected to 
complex thermomechanical loadings, such as those found in 
heat engine applications. CARES uses results from 
MSC/NASTRAN or ANSYS finite-element analysis programs 
to evaluate how inherent surface and/or volume type flaws 
affect component reliability. CARES utilizes the Batdorf model 
and the two-parameter Weibull cumulative distribution 
function to describe the effects of multiaxial stress states on 
material strength. The principle of independent action (PI A) 
and the Weibull normal stress averaging models are also 
included. Weibull material strength parameters, the Batdorf 
crack density coefficient, and other related statistical quantities 
are estimated from four-point bend bar or uniform uniaxial 
tensile specimen fracture strength data. Parameter estimation 
can be performed for a single or multiple failure modes by 
using a least-squares analysis or a maximum likelihood 
method. Kolmogorov-Smirnov and Anderson-Darling 
goodness-of-fit tests, 90-percent confidence intervals on the 
Weibull parameters, and Kanofsky-Srinivasan 90-percent 


confidence band values are also provided. Examples are 
provided to illustrate the various features of CARES. 

Introduction 

The unique properties that advanced ceramics offer in the 
areas of high-temperature strength, environmental resistance, 
and low density provide the potential for greatly increased fuel 
efficiency and reduced emissions in aerospace and automotive 
engine applications. Consequently, research has focused on 
improving ceramic material processing and properties as well 
as on establishing a sound design methodology. 

Because of the variable severity of inherent flaws, the nature 
of ceramic failure is probabilistic and optimization of design 
requires the ability to accurately determine a loaded com- 
ponent’s reliability. Methods of quantifying this reliability and 
the corresponding failure probability have been investigated 
and refined by the authors. The result of this effort is a public 
domain computer program with the acronym CARES 1 
(ceramics analysis and reliability evaluation of structures). 
CARES (ref. 1) calculates the fast fracture reliability of 
macroscopically isotropic ceramic components. These 

formerly SCARE (structural ceramics analysis and reliability evaluation). 
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components may be subjected to complex thermomechanical 
loadings such as those found in heat engine applications. 

The design methodology used by CARES combines three 
major elements: (1) linear elastic fracture mechanics (LEFM) 
theory which relates the strength of ceramics to the size, shape, 
and orientation of critical flaws, (2) extreme value statistics 
to obtain the characteristic flaw size distribution function, 
which is a material property, and (3) material microstructure. 
Inherent in this design procedure is that component integrity 
is a function of the entire field solution of the stresses and is 
not based only on the most highly stressed point. In addition, 
the size of the stressed material surface area and volume affect 
the component strength. 

Probabilistic component design requires the determination 
of the fracture strength distribution from simple geometry 
flexural or tensile test specimens. The statistical material 
parameters are estimated as functions of temperature, specimen 
loading, and geometry. From these data the reliability for a 
complex component geometry and loading is then predicted 
(fig. 1). Appropriate design changes are made until an 
acceptable probability of failure has been reached. 


•CERAMICS ARE BRITTLE AND HAVE MANY FLAWS 

•RANDOM FLAW SIZE AND ORIENTATION REQUIRE PROBABILISTIC METHOD 

• APPROACH: 



SIMPLE TESTS COMPLEX PREDICTIONS 


•REQUIRES ENTIRE STRESS FIELD, NOT MAXIMUM STRESS POINT 

Figure 1.— CARES uses statistical material strength parameters determined 
from simple tests to predict the reliability of complex ceramic components. 


Program Capability and Description 

CARES is an integrated computer program written in 
FORTRAN 77 which uses Weibull (refs. 2 and 3) and Batdorf 
(refs. 4 and 5) fracture statistics to predict the fast fracture 
reliability of isotropic ceramic components. CARES has three 
primary functions: (1) to analyze statistically the data obtained 
from the fracture of simple uniaxial tensile or flexural 
specimens, (2) to estimate the Weibull and Batdorf material 
parameters using these data, and (3) to perform a fast fracture 
reliability evaluation of a ceramic component experiencing 
thermomechanical loading. Component reliability is predicted 
using elastostatic finite-element analysis output from the 
MSC/NASTRAN or ANSYS computer programs. 

The CARES code includes a number of fracture theories 
to predict material response due to multiaxial stresses. These 
methods are summarized in table I. The Batdorf (ref. 5) method 
is recommended because it couples LEFM with the Weibull 
weakest link theory (WLT) (ref. 2). The Weibull normal stress 
averaging method (ref. 3) and the principle of independent 
action (PI A) (refs. 6 and 7) theories are included for 
comparison purposes and because of their previous popularity. 
All the fracture models available in CARES use the Weibull 
two-parameter probability of strength distribution. 

Figure 2 shows the fracture criteria and crack geometries 
available to the user for both surface and volume distributed 
flaws. The PI A and Weibull normal stress averaging fracture 
theories do not require a crack geometry. Batdorf s fracture 
theory can be used with several different mixed-mode fracture 
criteria and crack geometries. For coplanar crack extension, 
CARES uses the total strain energy release rate theory. Out- 
of-plane crack extension criteria are approximated by a simple 
semiempirical equation (refs. 8 and 9). This equation involves 
a parameter which can be used to approximate various mixed- 
mode theories or experimental results. For comparison, 
Griffith’s maximum tensile stress analysis for volume flaws 
is also included. The highlighted boxes in figure 2 show the 
recommended fracture criteria and flaw shapes. 

The statistical material parameters are obtained from the 
fracture stresses of many test specimens (ideally 30 or more) 
of fixed geometry and loading. Solutions for the four-point 


TABLE I.— STATISTICAL FAST FRACTURE MODELS AVAILABLE WITH CARES 


Weakest link 
fracture model 

Size 

effect 

Stress 

state 

effects 

Computa- 

tional 

simplicity 

Theoretical basis 

Weibull (1939) 

Yes 

Uniaxial 

Simple 

Phenomenological 

Normal stress 



Multiaxial 

Complex 

Phenomenological 

averaging (1939) 






Principle of 



Multiaxial 

Simple 

Maximum principal 

independent action (1967) 





stress theory 

Batdorf: 



Multiaxial 

Complex 

Linear elastic 

Shear-insensitive (1974) 





fracture mechanics 

Shear-sensitive (1978) 
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Figure 2.— Available failure criteria and crack shapes. (Recommended 
failure criteria and crack shapes are highlighted.) 


modulus of rupture (MOR) bending bar (ref. 10) and the pure 
tensile specimen (ref. 1 1) tested at a user specified temperature 
have been incorporated in the CARES program. Since the 
statistical material parameters are a function of temperature, 
up to 20 data sets may be input at discrete temperature levels. 
Lagrangian polynomials are utilized to interpolate the 


parameter values at other temperatures. Each data set may 
consist of up to 200 specimen fracture stresses. Each specimen 
can be identified by its mode of failure— either volume flaw, 
surface flaw, or unknown— so that statistical parameters from 
competing failure modes can be estimated. 

CARES can identify potential bad fracture stress data 
(outliers). The outlier test developed by Stefansky (ref. 12) 
and subsequently used by Neal et al. (ref. 13) is employed. 
The test can detect multiple outliers from a sample of up to 
100 specimens at either the 1-, 5-, or 10-percent significance 
levels. Data detected as outliers are flagged with a warning 
message and any further action is left to the discretion of the 
user. 

Weibull parameters are obtained as a function of the 
specimen surface and/or volume as requested by the user. 
Biased estimates of the Weibull parameters are obtained from 
either least-squares analysis or the maximum likelihood method 
for complete or censored samples (competing failure modes). 
CARES uses the Weibull log likelihood equations given in 
Nelson (ref. 14) and the rank increment adjustment method 
described by Johnson (ref. 15). 

Because the estimates of statistical parameters are obtained 
from a finite amount of data, they contain an inherent 
uncertainty which can be characterized by bounds in which 
the true parameters are likely to lie. For the maximum 
likelihood method the 5- and 95-percentile confidence limits 
for the Weibull parameters are provided (ref. 16). 

The ability of the probability distribution calculated from 
the Weibull parameters to reasonably fit the experimental 
data is measured with the Kolmogorov-Smirnov (K-S) and 
Anderson-Darling (A-D) goodness-of-fit tests. These tests are 
discussed by D’Agostino and Stephens (ref. 17). The A-D 
test is provided because it is more sensitive to discrepancies 
at low and high probabilities of failure than the more 
commonly used K-S test. The Kanofsky-Srinivasan 90-percent 
confidence band values (ref. 18) about the Weibull line are 
given as an additional test of the fit of the data to the Weibull 
distribution. 

Figure 3 illustrates the operational flow of the program. For 
a finite-element model, reliability calculations are performed 
at the element level and the overall component reliability is 
then the product of all the element survival probabilities. 
Reliability due to the presence of volume flaws is calculated 
from the volume statistical material strength parameters and 
the output of the stresses, volumes, and temperatures from 
isoparametric brick, wedge, or axisymmetric finite elements. 
Reliability evaluation due to the presence of surface flaws is 
calculated from the surface statistical material strength 
parameters and the output of the two-dimensional surface 
stresses, areas, and temperatures from isoparametric quad- 
rilateral and triangular shell elements. Solid elements are used 
for the structural modeling. Shell elements are only used to 
identify external surfaces of solid elements consistent with the 
component external boundaries that are required for the 
reliability analysis. 
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Figure 3.— Block diagram for analysis and reliability evaluation of ceramic 

components. 

Provision is made in CARES to permit the use of cyclic 
symmetry modeling. CARES also has a multiple material 
capability; a model can consist of up to 100 different materials 
(up to 100 different statistical material characterizations). 
Elements not designated as brittle materials are ignored in the 
reliability computations. Temperature-dependent statistical 
material properties are interpolated at each individual element 
temperature. The risk of rupture intensity is also calculated 
for each element and these values are sorted to determine the 
maximum values. 

Two versions of the code, designated as CARES 1 and 
CARES2, are available. The CARES1 version assumes that 
stress and temperature gradients within each element are 
negligible, and, therefore, only element centroidal principal 
stresses are used in the reliability calculations. The CARES2 
version takes into account element stress gradients by dividing 
each brick element into 27 subelements and each quadrilateral 
shell element into 9 subelements. Subelement centroidal 
principal stresses are then computed and used in the subsequent 
reliability calculations. CARES2 enables the finite-element 


model to consist of fewer elements for the same level of 
convergence to the true solution as CARES 1. 

Input Requirements 

To control the execution of the CARES program, an input 
file must be prepared. On the tape or disks provided with the 
program is a file called TEMPLET INP that can be used to 
construct an input file for a particular problem. Input to 
CARES is keyword driven. Data are input by the user under 
each keyword. An explanation of the input required or a list 
of input choices is provided next to the keyword. 

The CARES program requires three categories of input: (1) 
Master Control Input, (2) Material Control Input, which 
includes temperature-dependent material data, and, optionally, 
(3) MSC/NASTRAN or ANSYS output data files from finite- 
element analysis. The Master Control Input is a set of control 
indices that directs the overall program execution. The 
Material Control Input consists of control indices and either 
the data required to estimate the statistical material parameters 
or direct input of the statistical parameter values themselves 
for various temperatures. This input category includes the 
choices of fracture criteria and flaw shapes shown in figure 2. 
The Master Control Input and the Material Control Input are 
contained in the TEMPLET INP file. The third input category, 
MSC/NASTRAN or ANSYS output data files, includes finite- 
element analysis data files containing the element stresses, 
volumes/areas, and temperatures. 

Output Information 

The first part of the CARES output is an echo of the choices 
selected (or default values) from the Master Control Input. 
If a finite-element model reliability analysis is not performed, 
then CARES proceeds to echo the Material Control Input. If 
postprocessing of a finite-element model is done, then, for each 
element, the centroidal or subelement principal stresses with 
appropriate element area or volume and temperature are listed. 
The printing of element stress tables in CARES is optional. 
In addition, two element cross-reference tables are printed. 
The first table lists the shell element number and gives the 
corresponding solid element to which it is attached. The second 
table lists the solid element identification number and lists up 
to six associated shell elements (a brick element could have 
all of its six faces as external surfaces). 

CARES echoes the user inputs for each section of the 
Material Control Input. If statistical material parameters are 
directly input, then output pertaining to calculated values of 
the normalized Batdorf crack density coefficient k B will 
follow. If statistical material parameters are determined from 
experimental fracture data, then the output will identity the 
method of solution, the number of specimens in each batch, 
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and the temperature of each test. In addition, the output echoes 
the sorted input values of all specimen fracture stresses with 
proper failure mode identification. 

Results from the statistical analysis of the fracture data are 
then printed. The fracture strength and corresponding 
significance level are listed for detected outliers followed by 
the estimated statistical material parameters from least-squares 
or maximum likelihood analyses. The biased and the unbiased 
values of the Weibull shape parameter, the specimen Weibull 
characteristic strength, the upper and lower bound values at 
90-percent confidence level for both parameters, the specimen 
Weibull mean value, and corresponding standard deviation are 
printed for each specified temperature. For censored statistics 
these values are generated first for volume flaw analysis and 
subsequently for the surface flaw analysis. It should be noted 
that not all the previously mentioned information is available 
for all methods of material parameter estimation. 

The K-S goodness-of-fit test is done for each specimen 
fracture stress and the corresponding K-S statistics D+ , D— , 
and significance level are listed. Similarly, the K-S statistic 
for the overall sample set is printed along with the significance 
level. This overall statistic is the absolute maximum of 
individual data D+ and D— factors. For the A-D goodness- 
of-fit test, the A-D statistic A 2 is determined for the overall 
population and its associated significance level is printed. 

The next table of the output contains data to construct K-S 
90-percent confidence bands about the Weibull distribution. 
The table includes fracture stress data, the corresponding 
Weibull probability of failure values, the 90-percent upper and 
lower confidence band values about the Weibull line, and the 
median rank value for each data point. 

The last table from the statistical analysis section of CARES 
summarizes the material parameters used in component 
reliability calculations. These parameters, which are listed as 
a function of temperature, include the biased Weibull modulus, 
the normalized Batdorf crack density coefficient, and the 
Weibull scale parameter or the unit volume or unit area 
characteristic strength, whichever is appropriate. The values 
printed correspond to the experimental temperatures input and 
five additional interpolated sets of values between each input 
temperature. The interpolated parameters are output for 
checking purposes. Information on the selected fracture 
criterion and the crack shape is printed as required. 

If a component reliability analysis with finite-element data 
is being performed, then tables will be generated to summarize 
the reliability evaluation of each finite element. One table is 
provided for volume flaw analysis (solid elements) and one 
table is given for surface flaw analysis (shell elements) as 
requested by the user. The tables list the element identification 
(ID) numbers and the corresponding element material ID, 
survival probability, failure probability, risk-of-rupture 
intensity (risk-of-rupture divided by element volume or area), 
and temperature-interpolated statistical material parameters. 
Following each table is a sorted list of the 15 most critical 


risk-of-rupture intensity values and corresponding element 
numbers. Also included is the probability of failure and 
survival for the component surface or volume, whichever is 
appropriate. Finally, the overall component probability of 
failure and the component probability of survival are printed. 


Theory 

The first probabilistic approach used to account for the 
scatter in fracture strength of brittle materials was introduced 
by Weibull (refs. 2 and 3). This approach is based on the 
previously developed weakest link theory (WLT) (refs. 19 to 
21), which is primarily attributed to Pierce who proposed it 
while modeling yarn failure. The WLT is analogous to pulling 
a chain where catastrophic failure occurs when the weakest 
link in the chain breaks. The reliability of the chain is the 
product of the survival probabilities of the individual links. 

Phenomenological observations indicate that monolithic 
ceramic failures behave in accordance with WLT. For a 
ceramic component containing volume flaws and loaded in 
uniaxial tension, the probability of survival is expressed as 


P sV = exp 


\ v Ny{c)dv] 


( 1 ) 


where V is the component volume and the subscript V denotes 
volume-dependent terms. The function N v {o), referred to as 
the crack density function, represents the number of flaws per 
unit volume having a strength equal to or less than a. The value 
of the integral is called the risk-of-rupture. 

Weibull introduced a three-parameter power function for 
the crack density function N v {o)\ 


A I v {a) = 



( 2 ) 


where o uV is the threshold stress (location parameter), usually 
taken as zero for ceramics. The location parameter is the value 
of applied stress below which the failure probability is zero. 
When the location parameter is zero, the two-parameter 
Weibull model is obtained. The scale parameter o oV then 
corresponds to the stress level where 63.2 percent of specimens 
with unit volume would fracture. The scale parameter has 
dimensions of stress X (volume) 1/m L The reciprocal of o™y 
is called the uniaxial Weibull crack density coefficient k wV , 
and m v is the shape parameter (Weibull modulus or Weibull 
slope), which is a measure of the degree of strength dispersion. 
These statistical parameters are material properties that are 
temperature and processing dependent. They are evaluated 
from uniaxial flexure or tensile specimen fracture data. 

For surface flaw induced failure in ceramic structures these 
expressions become a function of the component surface area 
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A . Herein the subscript S denotes analogous parameters that 
are a function of surface area. 

To predict material response under multiaxial stress states, 
Weibull (ref. 3) proposed the normal stress averaging method. 
While this approach is intuitively plausible, it is somewhat 
arbitrary. Subsequently, Barnett and Freudenthal (refs. 6 and 
7) proposed the PI A model. The PI A fracture theory is the 
weakest link statistical equivalent of the maximum stress failure 
theory and is only applicable for tensile states of stress. The 
Weibull method of averaging the tensile normal stress and the 
PIA model have been the most widely used methods for brittle 
material design. However, they do not specify the nature of 
the defect causing the failure. 

Recognizing that brittle fracture is governed by LEFM, 
Batdorf (refs. 4 and 5) proposed that reliability predictions 
should be based on a combination of the weakest link theory 
and fracture mechanics. Conventional fracture mechanics 
dictates that both the size of the critical crack and its orientation 
relative to the applied loads determine the fracture stress. 
However, in brittle ceramics the small critical flaw size and 
the large number of flaws prevent determining the critical flaw, 
let alone determining its size and orientation. Instead, the 
combined probability of the critical flaw being within a certain 
size range and being oriented so that it may cause fracture 
is calculated. As flaw sizes correspond to strength levels and 
since strength is easier to measure than size for these 
microscopic flaws, the probability of a crack existing within 
a critical strength range is determined. 

The Batdorf theory assumes random flaw orientation and 
a consistent crack geometry. The component failure probability 
for volume flaws is expressed as 

P f v= 1- 


which depends only on the specimen fracture data, the Batdorf 
coefficient requires a fracture criterion and crack shape. 

If a shear-insensitive condition is assumed, fracture occurs 
when o n = o e > a cn where a n is the normal tensile stress on 
the flaw plane. However, for a flat crack it is known from 
fracture mechanics analysis that a shear stress r applied parallel 
to the crack plane (mode II or III) also contributes to fracture. 
Therefore, for polyaxial stress states, expressing the effective 
stress o e as a function of both a n and r is more accurate than 
assuming shear insensitivity. 

The equations derived by Batdorf and Heinisch (ref. 5) are 
based on self-similar (coplanar) crack extension. However, 
a flaw experiencing a multiaxial stress state usually undergoes 
crack propagation initiated at some angle to the flaw plane 
(noncoplanar crack growth). Shetty (ref. 9) performed experi- 
ments on polycrystalline ceramics and glass considering crack 
propagation as a function of an applied far-field multiaxial 
stress state. He modified an equation proposed by Palaniswamy 
and Knauss (ref. 8) to empirically fit experimental data. Using 
this criterion and a penny-shaped crack we obtain 



where v is Poisson’s ratio and C is a shear sensitivity constant 
adjusted to best fit the data. Shetty (ref. 9) found a range of 
values (0.80 < C < 2.0) for the materials he tested which 
contained large induced flaws. Note also that the CARES code 
has other fracture criteria and flaw shapes available, as 
indicated in figure 2. 

For mixed-mode fracture due to surface flaws the Batdorf 
failure probability equation is 


exp 


w max £2 (D,q cr ) dN v (a cr ) 
47 r do cr 

JkJo 


where cr cr , the critical stress, is defined as the remote, 
uniaxial, fracture strength of a given crack in mode I loading. 
The solid angle fi(E,cr cr ) is the area of a unit radius sphere 
containing all the crack orientations for which o e > o cr due 
to the existing stress state £. The effective stress a e is defined 
as the equivalent mode I stress on the flaw. The constant 47 t 
is the surface area of a unit radius sphere and corresponds to 
a solid angle containing all possible flaw orientations. The limit 
of integration o €m ^ is the maximum effective stress. The 
Batdorf crack density function N v (a cr ) is approximated by the 
power function 


N V {0 cr ) = k BV C 


(4) 


Pfs= 1 - 


exp 


e ma\ co(E,a cr ) dN s (o cr ) 


A JO 


27T 


do ri 


do cr dA 


( 6 ) 


where u(L,o cr ) is the total arc length on a unit radius circle 
on which the projection of the equivalent stress satisfies 
o e > o cn and 27 t is the perimeter of the circle. The cracks are 
assumed to be randomly oriented in the plane of the external 
boundary with their planes normal to the surface. 

For noncoplanar crack growth and a semicircular surface 
crack we obtain for the effective stress 



where the material Batdorf crack density coefficient k BV and 
Weibull modulus m v are evaluated from experimental uni- 
axial fracture data. In contrast to the Weibull coefficient k wV , 


where again C is adjusted to best fit the data. 

Selected statistical theories and equations for parameter 
evaluation are explained in detail in reference 22. Typically, 
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for brittle materials, the Weibull parameters are determined 
from simple specimen geometry and loading conditions, such 
as beams under flexure and either cylindrical or flat specimens 
under uniform uniaxial tension. The flexural test failure 
probability can be expressed in terms of the extreme fiber 
fracture stress oy or modulus of rupture MOR using the two- 
parameter Weibull form 



where m is the Weibull modulus and a e is the volume or area 
specimen characteristic strength (characteristic modulus of 
rupture, MOR 0 ). The Weibull scale parameter a 0 for volume 
and surface cracks is determined from o e , m , and the 
specimen geometry and loading. The characteristic strength 
a d is defined as the uniform stress or extreme fiber stress at 
which the probability of failure is 0.6321. 

Example 1 — Statistical Material 
Parameter Estimation 

To illustrate the methods used to estimate statistical material 
parameters, results from the fracture of four-point bend bars 
broken at NASA Lewis Research Center and analyzed by 
CARES are compared to results independently obtained by 
Bruckner-Foit and Munz (ref. 23) for the International Energy 
Agency (IE A) Annex II agreement (ref. 24). Two different 
materials were analyzed— namely, a hot isostatic pressed 
(HIPped) silicon carbide (SiC) from Elektroschmelzwerke 
Kempten (ESK), West Germany, and a HIPped silicon nitride 
(Si 3 N 4 ) from ASEA CERAMA, Sweden. 

The flexhre bars were distributed by Oak Ridge National 
Laboratory (ORNL) to the five participating U.S. laboratories, 
including NASA Lewis. The bars were fractured at these 
laboratories, and the fracture stress data sets were returned 
to ORNL as complete data without censoring for different 
failure modes. The number of specimens of a particular 
material given to each U.S. participant was 80. 


The results of the 80 silicon carbide flexure bars tested at 
NASA Lewis were analyzed using the CARES code to 
calculate the least squares and maximum likelihood estimates 
(MLE’s) of the Weibull parameters. The MLE values from 
CARES for a complete sample are compared in table II to the 
values obtained from reference 23. The SiC fracture data are 
plotted in figure 4 along with the proposed Weibull line and 
the Kanofsky-Srinivasan 90-percent confidence bands. Since 
no outliers were detected, all the data are within the 90-percent 
bands, and the goodness-of-fit significance levels are high, it 
is concluded that the fracture data show good Weibull 
behavior. 

For the ASEA CERAMA HIPped Si 3 N 4 bars fractured at 
NASA Lewis and subsequently analyzed as a complete sample, 
the statistical material parameters estimated with CARES and 
reference 23 are shown at the top of table III. The comparison 
of the MLE’s with reference 23 is very good. When analyzed 
by CARES as a complete sample, the significance levels of 
54 and 35 percent from the K-S and A-D goodness-of-fit tests, 
respectively, were relatively low, indicating a questionable fit 
to the proposed Weibull distribution. From the outlier test the 
highest strength fracture stress was detected to be an outlier 
at the 1 -percent significance level. Several of the lower 
strengths were flagged as outliers at various significance levels 
(1, 5, or 10 percent). Figure 5 shows a Weibull plot of the 
data. In this figure it appears that the data are bimodal with 
an outlier point at the highest strength. 

Because of the observed trends, the data were reanalyzed 
assuming a censored distribution and removing the highest 
strength value (oy= 817.2 MPa) as bad data. Although it is 
possible that both failure modes were surface induced, for the 
sake of this example it is assumed that the low strength failures 
were predominantly due to volume flaws and the high strength 
specimens predominantly fractured due to surface flaws. Since 
results from fractography of the individual specimens to 
identify the various failure modes were not available, the 
fracture origins had to be arbitrarily assigned prior to param- 
eter estimation. Note that identifying individual specimen flaw 
origins is especially important for small sample sizes where 
a plot of the data does not yield clear trends. However, for 


TABLE II. -WEIBULL PARAMETERS, CONFIDENCE LIMITS, KOLMOGOROV-SMIRNOV, AND ANDERSON-DARLING 
TEST RESULTS FOR SILICON CARBIDE FOUR-POINT BEND BAR FRACTURE DATA 

[All estimates are biased estimates; 80 samples per material; complete sample analysis.] 


Method 

of 

analysis 

Source 

of 

data 

Shape 

parameter, 

m 

90-Percent 
confidence 
limits on 
m 

Characteristic 

strength, 

MPa 

90-Percent 
confidence 
limits on 

°d> 

MPa 

K-S test 
statistic, 
D 

K-S test 
significance 
level, 
a, 

percent 

A-D test 
significance 
level, 
a, 

percent 

Upper 

Lower 

Upper 

Lower 

Maximum likelihood 

CARES 

6.48 

7.38 

5.52 

556 

573 

539 

0.070 

83 

86 

Maximum likelihood 

Ref. 23 

6.59 

7.65 

5.61 

556 

574 

539 

.063 

(a) 

(a) 

Least squares 

CARES 

6.59 

— 

— 

555 

— 

— 

.071 

81 

83 


a Not available. 
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Figure 4.— Weibull distribution and 90-percent confidence bands determined 
from fracture stress data for ESK HIPped silicon carbide four-point bend 
bars broken at NASA Lewis (not all data points shown). 



Figure 5.— Complete sample, assumed censored sample Weibull distributions, 
and outlier determined from fracture stress data for ASEA CERAMA 
HIPped silicon nitride four-point bend bars broken at NASA Lewis (not 
all data points shown). 


TABLE III. -WEIBULL PARAMETERS, KOLMOGOROV-SMIRNOV, AND ANDERSON-DARLING TEST 
RESULTS FOR SILICON NITRIDE FOUR-POINT BEND BAR FRACTURE DATA 


Method 

of 

analysis 

Source 

of 

data 

Assumed 

distribution 

Sample 

size 

Shape 

parameter, 

m 

Characteristic 

strength, 

MPa 

K-S test 
statistic, 
D 

K-S test 
significance 
level, 
a, 

percent 

A-D test 
significance 
level, 
a, 

percent 

Maximum likelihood 

CARES 

Surface 

80 

13.39 

686 

0.0901 

54 

35 

Maximum likelihood 

Ref. 23 



80 

13.40 

686 

.088 

(a) 

(a) 

Least squares 

CARES 



80 

11.74 

691 

.128 

15 

11 

Maximum likelihood 





79 

16.22 

683 

.078 

73 

56 

Least squares 





79 

11.98 

688 

.124 

18 

10 

Maximum likelihood 



Volume 

9 

4.13 

1128 

.081 

69 

58 




Surface 

70 

22.81 

692 

.081 

69 

58 

Least squares 



Volume 

9 

6.74 

830 

.112 

28 

35 




Surface 

70 

22.93 

691 

.112 

28 

35 

Maximum likelihood 



Volume 

13 

6.79 

876 | 

.074 

78 

88 




Surface 

66 

21.00 

693 

.074 

78 

88 

Least squares 



Volume 

13 

6.84 

864 

.085 

62 

38 




Surface 

66 

15.87 

697 

.085 ! 

62 

38 


a Not available. 


the NASA Lewis Si 3 N 4 data the sample size was large and 
clear trends could be observed, although extra care is required 
to determine if the trends are surface or volume flaw based. 

Two censored distributions were analyzed. The first was 
based on an inspection of figure 5, where the lowest nine 
strengths were assumed due to volume flaws and the remainder 
due to surface flaws. The second assumed thirteen volume 
flaws, where the particular volume flaw fracture strengths were 
assigned such that the MLE’s more closely fit the experimental 
data. 


Results of this analysis are shown in table III. At the top 
of the table are the parameter estimates and goodness-of-fit 
statistics for the 80-specimen complete sample followed by 
a complete and censored sample analysis of the data with the 
outlier removed. The table shows the MLE’s along with the 
least-squares estimates. 

When the goodness-of-fit scores are used as the basis of 
judgement, the censored sample estimated using the maximum 
likelihood analysis and thirteen volume flaws gives the best 
fit to the data. The high scores of 78 and 88 for the K-S and 
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A-D tests, respectively, indicate good Weibull behavior. 
Improvements in the goodness-of-fit scores may be gained 
through correctly identifying the location of fracture origins. 

It should be noted from figure 5 that the assumed volume 
flaw distribution dominates the failure response at low 
probabilities of failure. Therefore, in component design, it is 
essential to properly account for competing failure modes, 
otherwise nonconservative design predictions can result. 

Example 2— Rotating Annular Disk 

Failure predictions using the available fracture models from 
CARES, for both volume and surface analysis, were compared 
to failure predictions obtained by Swank and Williams (ref. 25) 
for a silicon nitride annular disk rotating at various speeds 
(fig. 6). The Weibull material parameters were independently 
evaluated from four-point MOR bar tests using a total of 85 
specimens. The Weibull modulus m was 7.65 and the 
characteristic modulus of rupture MOR 0 was 808 MPa. If 
fracture was assumed to be due to volume flaws, then 
a oV = 75 MPa (/n)° 3922 ; if fracture was assumed to be due to 
surface flaws, then, a oS = 232 MPa (w)° 2614 . Swank and 
Williams assumed that both the bars and disks broke because 
of volume flaws. 

Seven disks were fracture tested, and the experimental disk 
Weibull modulus of 4.95 was considerably different than the 
7.65 value based on MOR specimen data. A better agreement 
between disk and MOR Weibull slopes would lead to improved 
predictions in failure probabilities for the entire data range. 
Estimating parameters from small sample sets greatly increases 
potential deviation from the true population parameters. 
Confidence limits are used to measure the intrinsic uncer- 
tainties in parameter estimates from finite sample sizes. If 
potential experimental errors are excluded, it is possible that 
the rotating disks and the MOR bars broke due to the same 



Figure 6.— Rotating annular disk with 15° sector finite-element mesh con- 
taining eight brick elements (not to scale). Material, NC-132 hot-pressed 
Si 3 N 4 ; inner disk radius, r,, 6.35 mm; outer disk radius, r 0 , 41.28 mm; 
disk thickness, /, 3.80 mm; RPM range, 70 000 to 114 000. 


flaw population because their 90-percent confidence limits 
overlap between 6.56 and 6.98. 

It should be noted that Swank and Williams (ref. 25) also 
spin tested and performed a volume flaw reliability analysis 
on a contoured hub and a turbine blade ring geometry. Swank 
and Williams found that the Weibull moduli obtained from 
MOR bars cut from the hub and blade ring were in good 
agreement with the Weibull moduli obtained from the hub and 
blade ring spin tests, respectively. They also noted that the 
material parameters used for the annular disk reliability 
analysis were obtained under less tightly controlled conditions 
than those used with the other geometries. 

All reliability calculations were done with brick and 
quadrilateral shell elements. Because of symmetry, only 8 solid 
and 18 shell elements were used in one 15° sector for the model 
of the disk (fig. 6). Only one element spans both the thickness 
and circumferential directions. The shell elements were 
attached to solid elements consistent with the model external 
surfaces. 

Experimental results are plotted in figure 7 along with shear- 
insensitive (normal stress) and various shear-sensitive pre- 
dictions from CARES2 for the volume flaw analysis. Results 
were compared to data calculated by Swank and Williams, who 
used the Weibull normal stress averaging method and linear 
axisymmetric finite elements. The agreement between failure 
predictions is good, with the small discrepancy probably due 
to the different stress-volume data used in solving the reliability 
problem. 

When using the Batdorf model with CARES there are two 
methods available to the user to describe the material shear 
sensitivity. With the input parameter IKBAT set equal to zero, 
k B is calculated using only the o n stress component on the 
crack plane. This option assumes that mode I fracture is 
intrinsic to uniaxial loading. When IKBAT is set equal to one. 



Figure 7.— Comparison of experimental failure probabilities with those for 
various fracture models for a rotating annular disk (volume flaw analysis): 
m^ = 7.65; a oV = 15 MPa (m) 0 3922 . For IKBAT = 0, k BV = 16.30; for 
IKBAT = 1, k BV = 2.99 (only for Griffith crack, Shetty criterion, C = 0.82); 
G t , total strain energy release rate criterion. 
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then k B is calculated using the user selected fracture criterion 
and crack geometry. These two methods will yield opposite 
trends relative to the shear-insensitive criterion as shown in 
figure 7. When I KB AT is set equal to zero, the subsequent 
reliability predictions are more conservative. The value of 
I KB AT is chosen so as to best fit the reliability predictions to 
the experimental data. 

For this example the laboratory measurements agree best 
with the shear-sensitive fracture models using the I KB AT = 0 
option. Note that results are given for an approximation of the 
maximum strain energy release rate criterion G ^ using a 
Griffith crack, C = 0.82, and IKBAT = 1. Failure probabilities 
calculated from decreasing shear-sensitive effective stress 
equations move the probability of failure curves toward the 
shear-insensitive case. It is observed that Shetty’s criterion for 
C = 0.80 and IKBAT = 0 with the penny-shape crack gives 
the best agreement with experimental data. 

The disks were reanalyzed assuming that fractures in the 
MOR bars as well as the disks were caused by surface flaws. 
Selected results are shown in figure 8. The same trends are 
observed for the shear sensitive results as the volume analysis 
results. However, for a given speed, failure probabilities are 
significantly less than those obtained by the volume flaw 
analysis for all fracture models, indicating that material failure 
was most likely due to volume flaws. The main reason for 
the decreased failure estimates is the much higher equivalent 
surface Weibull scale parameter a oS . The importance of post- 
mortem fractography to identify the nature of the fracture 



Figure 8.— Comparison of experimental failure probabilities with those for 
various fracture models for a rotating annular disk (surface flaw analysis): 
m s = 7.65; a_ oS = 232 MPa (m) 026 ' 4 . For IKBAT = 0, k BS = 4.99\ for 
IKBAT =\,k BS =\.16 (only for Griffith crack, Shetty criterion, C = 0.82); 
G T , total strain energy release rate criterion. 


causing flaws is evident from the two different sets of answers 
in figures 7 and 8. 


Example 3 — Si3N4 Mixed Flow Rotor 

The CARES code is used by numerous companies world- 
wide in the automotive, aerospace, nuclear, and computer 



900 

Figure 9.— Silicon nitride mixed-flow rotor temperature distribution. 
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2 x 10-2 



PRINCIPLE OF INDEPENDENT 
ACTION, P f = 0.00082 


STRAIN ENERGY RELEASE RATE 
CRITERION, PENNY-SHAPED 
CRACK, P f = 0.0062 


1x10-2 -I 
1x10-4 

1x10-6 -I 



1 x io-io - I 


1 x 10-12 
1 x 10-14 



Figure 10.— Volume flaw analysis comparison of risk-of-rupture intensities for PIA (left) and total strain energy release rate criterion with a penny-shaped 
crack geometry and I KB AT = 0 (right). 


software fields. At NASA Lewis, CARES has been used for 
the preliminary design of a silicon nitride mixed-flow rotor 
for application in small, high-temperature engines. A single 
blade and a section of the rotor hub were analyzed using the 
cyclic symmetry option of MSC/NASTRAN. 

The results from the heat-transfer analyses are shown in 
figure 9. In figure 10 the element risk of rupture intensities 
are plotted from the CARES volume flaw analysis for the PIA 
and the total strain energy release rate criterion using 
IKBAT = 0 and the penny-shaped crack geometry. Note that 
the risk of rupture intensity is independent of the individual 
element geometry, unlike the probability of failure, and 
provides the designer with a means to visualize the critical 
stressed regions. 

Similar to the rotating annular disk example, the shear- 
sensitive criterion yields a higher probability of failure for the 
same applied load. However, the regions of low reliability are 
the same for both models. This design was optimized to yield 
a low probability of failure. It is observed in figure 10 that 
the most critical region is at the rotor hub. 


Conclusions 

The potential use of structural ceramics for high-temperature 
applications depends on the strength, toughness, and reliability 
of these materials. Components using ceramics can be designed 
for high reliability in service if the contributing factors that 
cause material failure are accounted for. This design 
methodology must combine the statistical nature of strength 
controlling flaws with fracture mechanics to allow for 
multiaxial stress states and concurrent flaw populations. This 
is accomplished using the NASA/CARES public domain 
computer program for predicting the fast fracture reliability 
of structural ceramic components. This framework will be 
subsequently built on to include ceramic fatigue due to 
subcritical crack growth. The CARES code and accompanying 
documentation may be obtained by contacting the authors. A 
PC-based version is also available for statistical analysis and 
parameter estimation only. 
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